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Resonant electromagnetic properties of nanoparticles fabricated from high-index semiconductor or 
dielectric materials are very promising for the realization of novel nanoantennas and metamaterials. In this 
paper we study optical resonances of Si nanocylinders located on a silica substrate. Multipole analysis of the 
experimental scattering spectra, based on the decomposed discrete dipole approximation, confirms 
resonant excitation of electric and magnetic dipole modes in the Si nanocylinders. Influences of light 
polarization and incident angle on the scattering properties of the nanocylinders are studied. It is shown that 
the dependence of resonant excitation of the electric and magnetic modes in the nanocylinders on incident 
angle and polarization of light allows controlling and manipulating the scattered light in this system. The 
demonstrated properties of Si nanocylinders can be used for the realization of dielectric metasurfaces with 
different functional optical properties. 



It is well-known from Mie theory that the first and second lowest resonances of dielectric particles correspond 
to the magnetic and electric dipole terms \ For higher refractive indexes, the quality factors of these Mie 
resonances and their light scattering efficiency increase. Due to these properties, the high-refractive index 
dielectric particles are considered as promising building blocks for the creation of metamaterials with negative 
effective permeability and permittivity^. Experimental demonstrations of dielectric metamaterials with strong 
magnetic and electric responses at microwave and mid-infrared frequencies were reported^'^. Recently, strong 
magnetic and electric dipole resonant responses of silicon spherical nanoparticles were predicted theoretically in 
the visible spectral range^ and in the infrared region^'^, and have been demonstrated experimentally^"^. Using 
these optical properties of Si nanoparticles new all-dielectric optical nanoantennas and optical metamater- 
^^jg5,9,i3,i4 ]^2ive been suggested. Recently, it has been theoretically^^ and experimentally^^ shown that the spectral 
positions of electric and magnetic dipole resonances of Si nanoparticles can be tuned by changing their shape. For 
example, it is possible to excite electric and magnetic dipole resonances at the same optical frequency for Si 
nanodisks by changing their aspect ratio^^'^^. This possibility allows controlling light scattering charcteristics by 
using Si nanodisks or nanocylinders with definite size parameters. This is one of the reasons why optical 
properties of cylindrical-shaped dielectric nanoparticles has been intensively studied recently^ It was experi- 
mentally demonstrated that a two-dimensional periodic array of silicon nanocylinders fabricated on a substrate 
and supporting Mie resonances yields almost zero total reflectance over the all visible spectral range^^. Optical 
resonances of cylindrical Si nanoparticles located on different substrates and irradiated by normally incident light 
waves have been considered theoretically^^. The influence of Mie resonances on the light absorption in periodic Si 
nanocylinders arrays has been experimentally demonstrated and studied in Ref. 20. In this paper, we investigate 
experimentally and theoretically optical resonant properties of single cylindrical Si nanoparticles located on a 
silica substrate and discuss the role played by excited multipole modes in their extinction and scattering spectra. 
Influences of light polarization and incident angle on the scattering diagram are studied as well. 

Results 

Experiment. Scanning electron microscope (SEM) images of fabricated Si nanocylinders located on a silica 
substrate are shown in Fig. 1. The height of all nanocylinders is determined by the thickness of Si layer on a 
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Figure 1 | (a) Dark-field images of cylindrical Si nanocylinders with 
diameters of — 109, —143, —203 nm (from the left to the right) on a silica 
substrate, (b) Corresponding top view SEM images of the nanoparticles 
and substrate. To improve contrast in SEM images the nanocylinders were 
covered by a thin gold layer ( — 15 nm) subsequent to dark- field 
measurements. The scale bar presents 100 nm and refers to all images in 
(b). 

Soitec SOI wafer (Fig. 9a) and is equal to 220 nm. Their diameters 
vary in the range of 109 to 203 nm. In order to increase the contrast 
of SEM images. Si nanocylinders were covered by a thin gold layer 
after measurements of their optical spectra. Corresponding SEM 
images of the Si nanocylinders with the gold covering are shown in 
Fig. lb. Dark-field images of the Si nanocylinders are presented in 
Fig. la. Different colors are produced by nanocylinders with different 
diameters indicating a strong dependence of scattering resonances 
on this parameter. Note that for nanocylinders with small diameters 
(109 nm and 143 nm) the dark- field images are homogeneous in 
color, whereas for the nanocylinder with diameter 203 nm the 
dark-field image has different colors. As it will be shown below, in 
the first case the dark- field images are determined by the resonantly 
excited dipole modes, in the second case the image is formed 
by several resonantly excited high-order modes at different 
wavelengths. 

In order to consider this effect in more details, we performed single 
nanoparticle spectroscopy measurements and corresponding 
numerical simulations. Figure 2 demonstrates the experimental 
back- scattering spectra of Si nanocylinders with different diameters. 
It is clear that all scattering spectra have several resonant peaks, 
which spectral positions depend on the particle diameters. The main 
(maximal) resonant peaks are red shifted with the increasing particle 
diameter. However, for the largest particle the main peak is in the 
near infrared and hence the particle visually appears due to higher 
order modes (Fig. la). The width and magnitude of the resonant 
peaks are also increasing with the growth of the diameters. Note that 
for the relatively big nanocylinder (diameter of 203 nm) in Fig. 2a 
there are several resonant peaks, whereas for nanocylinders with 
small diameters there is only single resonant peak (Fig. 2b). 

Calculated spectra. It is known that Si cylindrical nanoparticles can 
support optical resonances of both electric and magnetic 
^ggi5,i7,i8,2i Spectral positions of the resonances are dependent on 
the particle aspect ratio and incident-light direction and polarization. 
In order to clarify the role of Si- nanocylinder magnetic and electric 
modes in the resonances presented in Fig. 2 we calculated scattering 
cross-sections for Si nanocylinders using the method of decomposed 
discrete dipole approximation (DDDA). In our simulations we 
consider a single Si nanocylinder irradiated by linear-polarized 
plane wave which obliquely falls on a substrate with refractive 
index of 1.5 (Fig. 3), incident angle is chosen to be 70 degrees. 
Because in the experiment the incident light illuminates the sample 
from all sides the light polarization condition is not totally defined. 
Therefore, we shall consider interaction of scatterers with incident 
light waves of TM (the magnetic field of incident light is 
perpendicular to the XZ- plane in Fig. 3) and TE (the electric field 
of incident light is perpendicular to the XZ-plane in Fig. 3) 
polarizations. Dielectric permittivity for crystalline Si is taken from 
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Figure 2 | Experimental back- scattering spectra of cylindrical Si 
nanoparticles with diameters D. The particle heights are 220 nm. 

Ref 22. Optical spectra of scattering cross-sections calculated for Si 
nanocylinders with different diameters are presented in Fig. 4. 
Independent on light polarization, the resonant peaks are shifted 
to the red side with increasing nanocylinder diameters (Fig. 4a, b). 
The spectral behavior of the scattered cross -sections in Fig. 4 agrees 



Scattering angle 




Figure 3 | Schematic representation of the physical system involved in 
the numerical simulations: a cylindrical particle (height H and diameter 
D) placed on a glass substrate is illuminated by a light wave, only scattered 
waves propagating into the scattering angle-region give contributions 
into calculated scattering cross- sections. 
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with experimental curves presented in Fig. 2 indicating the adequacy 
of the theoretical model. The spectral positions of resonances are also 
independent on light polarization. For example, for the nanocylinder 
with D = 150 nm the main difference between the first resonances 
shown in Figs. 4a and 4b is only in relative values of its two shoulders 
(in the TE polarization case the second shoulder at >^ = 675 nm is 
more noticeable than in the case TM polarization). Such tendency is 
conserved for the nanocylinders with smaller diameters. However, 
for the nanocylinder with diameter D = 90 nm this difference is 
considerably reduced (Fig. 4c) and the scattering cross-sections 
agree with experimental spectra presented in Fig. 2b. In spite of 
qualitative agreement between the experimental spectra in Fig. 2 
and theoretical data in Figs. 4, there are noticeable differences. The 
diameters of nanocylinders in our experiment are larger than that in 
simulations, due to a small oxide shell created on Si nanoparticles 
during the fabrication process. Thus the effective size of the Si 
nanoparticle, which determines the resonance position in the 
experiment, decreases^. 

In order to clarify the origins of the main resonant peaks in Fig. 4 
we applied the multipole decomposition^^'^^ to the scattering electric 
fields and compared the scattering cross-sections, calculated sepa- 
rately for different multipole moments of a nanocylinder, with gen- 
eral scattering cross-section. In our simulations we calculated the 
multipole moments of the nanocylinder with respect to its center 
of mass. Results of this comparison are shown in Fig. 5 for the Si 
nanocylinder with diameter of 100 nm. The first low- frequency res- 
onance (wavelength = 590 nm) corresponds to the excitation of 
the magnetic dipole (MD) moment of the nanocylinder indepen- 
dently of light polarization. The shoulder in Fig. 5a and the peak 
Fig. 5b at same X = 540 nm are basically determined by the electric 
dipole (ED) moment of the nanocylinder. In the case of T£-polar- 
ization the MD scattering cross -section has two resonant peaks cor- 
responding to resonant excitation of the in-plane component {X = 
590 nm) and out-of-plane component {X = 520 nm) of the nanocy- 
linder's MD. Note that a magnetic quadrupole (MQ) resonance 
appears on the background of broad ED and MD resonances for 
both polarizations. Its influence on the calculated scattering spectra 
is relatively weak. Differences between positions of the resonant 
peaks of the total scattering cross-sections and the multipole scatter- 
ing cross- sections in Fig. 5 are explained by contributions of the 
interference between multipole fields into the total scattering 
cross-sections calculated for the scattering angle 80° as shown in 
Fig. 3. 

Regardless of a certain similarity between the scattering cross- 
sections for different light polarizations presented in Fig. 4 there 
are important differences with respect to multipole excitations. 
Under inclined light incidence at a large incident angle the main 
dipole components excited by light are out-of-plane (z-) component 
of ED and in-plane (x-) component of MD for the TM polarization 
and vice-versa for the TE polarization. Because the out-of-plane 
dipoles do not radiate light normally to substrate surface, the scatter- 
ing cross-section into the relatively small scattering angle (Fig. 3 and 
Fig. 4) is basically determined by the in-plane components of ED and 
MD. As a result we see similarities between corresponding scattering 
cross-sections presented in Figs. 4a and 4b. This situation is signifi- 
cantly changed if total extinction and full (into all directions) scatter- 
ing cross-sections are considered. From Fig. 6 one can see that the 
values of the extinction and scattering cross-sections in the T£-polar- 
ization case are much smaller than for the TM-polarization. In the 
latter case the main (out-of-plane) component of ED moment and 
the only (in-plane) component of MD moment are resonantly 
excited by light at the same wavelength region increasing the extinc- 
tion and scattering (Fig. 6a). In the case of TE polarization the reso- 
nances of the main (out-of-plane) component (1 = 520 nm) and 
in-plane component (1 = 590 nm) of MD moment resonance are 
reaUzed outside of the ED resonance (Fig. 6b). Moreover the con- 




400 



500 



600 

X (nm) 



700 



800 



0.018 
0.016 
0.014 
0.012 
0.01 
S 0.008 
0.006 
0.004 
0.002 



-B— D=150 nm 
-A— D=125 nm 
-a— D=100 nm 
D=90 nm 



TE (b) 




600 

X (nm) 



X 10 




400 450 500 



550 600 650 

X (nm) 



700 750 800 



Figure 4 | Calculated spectra of scattering cross- sections into the conical 
region with scattering angle 80 degrees (Fig. 3) of Si nanocylinders. (a) 

Scattering of TM-polarized light plane waves, (b) Scattering of TE- 
polarized light plane waves, (c) Scattering cross sections (SCS) of Si 
nanocylinder with diameter D = 90 nm calculated for TM- and TE- 
polarized incident light, the green curve was calculated for polarization 
averaging. The height of the all nanocylinders is equal to 220 nm, the 
incident angle is equal to 70 degrees. /I is the wavelength of the incident 
light. 
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Figure 5 | Calculated spectra of total scattering cross-sections into the 
conical region with scattering angle 80 degrees (Fig. 3) of a Si 
nanocylinder irradiated by light waves of (a) TM polarization and (b) TE 
polarization. The nanocylinder diameter is equal to 100 nm. Other 
parameters are presented in Fig. 4. The scattering cross sections calculated 
separately for different multipole moments of nanocylinder are also 
presented: ED - electric dipole; MD - magnetic dipole; EQ - electric 
quadrupole; MQ - magnetic quadrupole. 
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Figure 6 | Calculated spectra of total extinction and scattering cross- 
sections of a Si nanocylinder irradiated by light waves of (a) TM 
polarization and (b) TE polarization. The nanocylinder diameter is equal 
to 100 nm. Other parameters are presented in Fig. 4. Contributions of the 
nanocylinder electric dipole (ED) and magnetic dipole (MD) moments 
into the extinction cross-sections are also presented. 



tributions of the in -plane ED component into the extinction and 
scattering cross-sections in Fig. 6b are much smaller than the con- 
tributions of the out-of-plane ED component in Fig. 6a due to differ- 
ences between corresponding components of the ED polarizability 
tensor of cylindrical nanoparticles. Comparing the extinction and 
scattering cross -sections in Fig. 6 one can see that the relative light 
absorption is weaker for TM polarization. 

Influence of the incident angle. The extinction cross -sections of the 
Si nanocylinder in Fig. 5 calculated as a function of the incident angle 
and the wavelength are demonstrated in Fig. 7 for the two polari- 
zations. It can be seen that for small incident angles (from 0 to 30 
degrees) the extinction spectra are very similar for both polarizations. 
The main resonances excited by light with wavelengths of 520 nm 
and 590 nm correspond to the inplane electric and in-plane mag- 
netic dipole moments of the nanocylinder, respectively. For incident 
angles larger than 30 degrees there are principal differences. For TM- 
polarized light the cross-section at = 590 nm starts to increase 
achieving the maximum around the incident angle of 60 degrees 
(Fig. 7a). This effect is explained by the increase of both out-of- 
plane electric dipole and in-plane magnetic dipole moments of the 
scatterer. Note that the resonant wavelengths of these two dipole 



modes coincide, but the out-of-plane electric dipole resonance is 
significantly broader than the magnetic dipole resonance (Fig. 6a). 
For incident angles larger than 75 degrees all resonances become 
weak due to strong decreasing of total incident electric and 
magnetic fields (in this angle region the reflection coefficients 
approach to — 1 for both polarizations, therefore the total external 
electric and magnetic fields approach zero)^^. In contrast, in the case 
of TE polarization (Fig. 6b) the resonant values of the extinction 
(scattering) cross-sections decrease with increasing incident angle. 

Dependence of resonant excitation of the electric and magnetic 
dipole modes in the nanocylinders on incident angle and polariza- 
tion of light allows controlling and shaping the scattered light in the 
system. Figure 8 demonstrates the possibility to change the light 
scattering pattern by changing polarization under a fixed incident 
angle. For TM polarization there is a strong asymmetric scattering 
into the substrate (Fig. 8a) due to interference between the fields 
generated (scattered) by the resonantly excited out-of-plane electric 
and in-plane magnetic dipole moments of the Si nanocylinder. It 
presents Kerker's type scattering^^"^^ for Si nanocylinders located 
on a glass substrate. In the conditions of TM-polarization, for large 
incident angles, we have a situation that is qualitatively similar to the 
case when a Si nanodisk is irradiated normally by linear polarized 
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Figure 7 | Extinction cross-sections (/zm^) of Si nanocylinders from Fig. 5 
(height is 220 nm, diameter is 100 nm) calculated as a function of the 
incident angle and the wavelength for (a) TM-polarized and (b) TE- 
polarized incident light. Arrows in (a) show the dipole modes which give 
the main contributions into the extinction cross-sections. Dashed Unes 
correspond to the extinction spectra presented in Fig. 6. Arrows indicate 
the resonant contributions of corresponding muhipole moments. 



hght^^. In these conditions the electric and magnetic dipole reso- 
nances can be excited at the same spectral region (like it is shown 
in Fig. 6a) resulting in the highly directive scattering pattern^^ due to 
interference between the electric -dipole and magnetic -dipole 
scattered fields^^. Note that similar asymmetrical distrubution of 
electromagnetic field can be observed in the system of two metal 
nanoparticles with detuned electric-dipole responses^^. Symmetri- 
cal scattering pattern in case of TE polarization (Fig. 8b) is deter- 
mined by the in-plane component of the magnetic dipole moment. 
Black regions outside the white dash circles in Fig. 8, for 0 > nil, are 
the forbidden light zones^^ which are determined by the critical angle 
Oc = n — arcsin(l/n2) ~ 138°, where ^2 is the reflective index of the 
substrate. 

Discussion 

In summary, it has been experimentally demonstrated that cylin- 
drical Si nanoparticles with diameters in the range of 109 to 
203 nm and height of 220 nm positioned on a sihca substrate have 
resonant optical responses in the visible spectral range. Theoretical 
analysis, based on the decomposed discrete dipole approximation, 
has shown that the optical properties of the Si nanocylinders are 
determined by their magnetic and electric dipole modes resonantly 
excited by incident light. With increase of the nanocylinder's dia- 
meter the dipole resonances are shifted into the infrared region. The 
influence of propagation direction and incident light wavelength on 
the resonant excitation of electric and magnetic modes of Si nano- 
cylinders was considered. The main attention has been given to the 
case when the the electric and magnetic dipole modes are resonantly 
excited in the visible spectral range. It has been demonstrated that the 
total extinction and scattering cross-sections and scattering directiv- 
ity are strongly dependent on the orientation of Si nanocylinders 
with respect to the incident direction of light and its polarization. 
It was found that the directivity of scattered light can be manipulated 
by the polarization of incident light due to simultaneous resonant 
excitation of electric and magnetic dipole modes of Si nanocylinders. 
The demonstrated optical properties of Si nanocylinders can be used 
for the realization of multi-functional dielectric metasurfaces and 
optical devices for nanoantenna and sensing applications. Further 
investigations will be concentrated on optical properties of Si nano- 
cylinders based on resonant excitation of their higher-order multi- 
pole modes. 

Methods 

Experimental methods. The silicon cylinders were fabricated using an etching mask 
of ZEP-520 electron beam resist patterned by electron beam lithography and 
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Figure 8 | Scattering directivity (a.u.) presented in spherical coordinates calculated for the system presented in Fig. 3: a Si nanocylinder with 
height of 220 nm and diameter of 100 nm is irradiated by a light plane wave with wavelength of 590 nm, the incident angle is 65 degrees. The 

polar angle 0 e [0, 7i] is represented by the radial distance, whereas the azimuthal angle (p G [0; 2n] is displayed in the figure plane. The center of the circle 
6 = 0 corresponds to the normal scattering above substrate; the boundary of the circle 6 = n corresponds to the normal scattering into substrate. 
Scattering above (below) the substrate corresponds to the region inside (outside) the white dash circles corresponding to 6 = n/2. (a) TM-polarized 
incident light, (b) TE-polarized incident Ught. 
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Figure 9 | (a) Schematic presentation of the fabrication stages: initial stage shows a SOI sample before the lithography procedure; final stage shows a Si 
nanocylinder located on Si02 substrate, (b) Schematic presentation of the setup for spectroscopic measurements. 



subsequent inductively coupled plasma etching on a Soitec SOI wafer with a 220 nm 
thick Si layer on 2 /im buried oxide. The fabrication is schematically shown in Fig. 9a. 

Scattering spectra of the fabricated Si nanocylinders were measured using optical 
dark- field spectroscopy. The spectroscopic set-up (Fig 9b) included a BX51 micro- 
scope (Olympus) equipped with a halogen light source and a fibre-coupled grating 
spectrometer QE65000 (Ocean Optics). In these measurements a single Si nanopar- 
ticle is irradiated at highly oblique angles around 70 deg. and back scattered light 
collected from angles up to 64 deg. determined by the numerical aperture (NA) of the 
xlOO microscope objective MPLFLNBDP (Olympus). The microscope images were 
captured with a XC30 digital color camera (Olympus). 

Theoretical and numerical methods. For theoretical analysis, we calculated the 
extinction and scattering cross-sections of cylindrical Si nanoparticles located on a 
silica substrate using the method of decomposed discrete dipole approximation 
(DDDA). This method combines the ordinary discrete dipole approximation^''"^^ 
with calculations of cartesian multipole moments of electric dipole systems and 
their contributions into the extinction and scattering cross-sections. Details of this 
method can be found elsewhere^^'^\ Briefly, the scattering object is replaced by a cubic 
lattice of electric point dipoles with a polarizability dp. The corresponding dipole 
moment p^ induced in each lattice point j (with the radius-vector r^) is found by 
solving coupled- dipole equations. After solution of these equations, the induced 
polarization P(r) of the scatterer can be written as 



P(r) ^p(5(r-ro) - - QV(5(r-ro) + - [V x m(5(r-ro)] 

6 CD 

+ \ 6(VV^(r-ro)) - [V X MV(5(r-ro)] . 
6 2<ji 



(2) 



Here the polarization P is presented as a sum of multipole terms corresponding to 
point multipole moments (electric dipole moment p, electric quadrupole moment 
Q(ro), magnetic dipole moment m(ro), magnetic quadrupole moment M(ro), and 
electric octupole moment 6(ro)) located at a point Tq. It is convenient to choose the 
point To at the scatterer's center of mass^l Here we take into account only several first 
multipole terms considering small scatterers compared to the wavelength. The point 
multipole moments of the scattering object are simply calculated from the space 
distribution of ^^-^^ 



p = ^P,-; Q(ro)= ^Q^'Cro); m(ro) = ^ m,-(ro); 

j=l ;=1 ;=1 

N N 

M(ro) = 5] M(ro); 6(ro) = ^ ^■(ro), 

j=i j=i 



(3) 



Here (5(r — r^) is the Dirac delta function. Application of the multipole 
decomposition^^ to the polarization (1) allows to write 



(1) 



where the corresponding point multipole moments, associated with the single electric 
dipole pj, are under the sum symbols. Analytical expressions for these multipole 
moments are presented in Ref 21. Finally, the multipole decomposition (2) is used for 
calculations of the extinction and scattering cross-sections, as it is shown in 
Ref 21. 
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